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ABSTRACT 



We aim at understanding the massive star formation (MSF) limit m(r) = 
SlOM^r/pc) 133 in the mass-size space of molecular structures recently proposed by 
iKauffmann & Pillai I ( f2010h . As a first step, we build on the property that power-law density 
profiles for molecular clumps combined with a volume density threshold for overall star for- 
mation naturally leads to mass-radius relations for molecular clumps containing given masses 
of star-forming gas. Specifically, we show that the mass m c i um „ and radius r c i ump of molecular 
clumps whose density profile slope is p and which contain a mass m,/, of gas denser than 

a density threshold p,/, obeys: m c i ump = m p J^ ( A ^ p ' p h ^ ^clwnp- ^ n a seconc l ste P> we use 

the relation between the mass of embedded-clusters and the mass of their most-massive star 
to estimate the mimimum mass of star-forming gas needed to form a 1OM star. Assuming 
a star formation efficiency of SFE ~ 0.30, this gives m t h C rit — 150M©. In a third step, we 
demonstrate that, for sensible choices of the clump density index (p ~ 1 .7) and of the clus- 
ter formation density threshold («,/, ~ 10 4 cm 3 ), the line of constant m t /, cr j t ~ 150M ra in the 
mass-radius space of molecular structures equates with the MSF limit for spatial scales larger 
than 0.3 pc. Hence, the observationally inferred MSF limit of Kauffmann & Pillai is consistent 
with a threshold in star-forming gas mass beyond which the star-forming gas reservoir is large 
enough to allow the formation of massive stars. For radii smaller than 0.3 pc, the MSF limit is 
shown to be consistent with the formation of a 1OM star (mthcrit — 30M Q with SFE ~ 0.3) 
out of its individual pre-stellar core of density threshold n tn ~ 10 s cm -3 . The inferred density 
thresholds for the formation of star clusters and individual stars within star clusters match 
those previously suggested in the literature. 
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1 INTRODUCTION 

Massive stars play a crucial role in the evolution of galaxies. Their 
ionising radiations and stellar winds carve bubbles into their sur- 
rounding interstellar medium (I SM), thereby triggering the forma- 
tion o f new stellar generations dDale et alj|2009l : Deharveng et al.l 
l2O10l) . while their explosions a s supernovae disperse into the ISM 
their nucleosynthesis products dHeger et al. 112003b . Massive form- 
ing stars also provide us with the only signatures of star formation 
in other gala xies. They are th us important tracers of galaxy star for- 
mation rates dCalzetti 120081) . The local environment of our Galaxy 
is an obvious first place where we can secure a firm grasp on their 
formation conditions, that is, where do they form or not, and why? 
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On top of that, the mass of the most-massive star in star clusters 
and the star cluster mass (formally, the cluster stellar mass when 
it still resides in the molecular clump out of which it forms) ar e 
observed to be related to each other dWeidner & Kroupa II2OO0I) . 
The very nature of that relation remains strongly debated, how- 
ever. Is it a physical one, i.e. is there a cluster- mass-dependent limit 
upon the maximum stellar mass in a cluster dWeidner & Kroupal 
1 20061) ? Or, is it a statistical relation, i.e. does it merely repre- 
sent the random-sampling-driven average relationship between the 
masses of clusters and their most common maximum stellar masses 
dParker & Goodwin I [2007b ? In spite of this uncertainty, there is 
a wide consensu s that the vast majority of massive stars form 
in star clusters l|Lada & Lada 1 120031 ; IParker & Goodwin! 120071; 
IWeidner & Kroupa I l2006h~ The relation between the embedded- 
cluster mass, m ec i, and the mass of cluster's most-massive star, 
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m *,max, demonstrates that the mass of molecular gas forming a 
star cluster influences the mass of this cluster 's most-massive star 
since: 

= f(m ec i) = f(SFE x mcFRg) , (1) 

where mcFRg is the gas mass of the cluster-forming region (CFRg) 
at the onset of star formation and SFE its star formation efficiency. 
Insights in what determines the mass m CFRg of star-forming gas 
inside individual molecular clumps hold therefore the potential of 
shedding light on the mass of their respective most-massive star 
and, eventually, on how massive stars form. This will also help to 
model better the impact of stellar feedback upon molecular clumps 
and molecular clump luminosities, since both are sensitively stellar- 

mas s-dependent. 

iKauffmann et alj 1201 Obi ) recently performed a significant 
quantitative step towards a better understanding of massive star 
formati on conditions. Applying the 'dendrogram-technique' devel- 
oped bv lRosolowski et al] J2008h to column density observations of 
several solar neighbourhood molecular clouds, they measure the ef- 
fective radiu r, of many column density contours around column 
density peaks, as well as the projected mass, m(r), contained within 
each contour. In that manner, they obtain a sequence of mass- 
size measurements over a comprehensive range of spatial scales 
for eac h molecular cloud analysed (see fig. 1 in IKauffmann et al .1 
2010a. for an illustration of the method). Their combined analy- 
sis of the Taurus, Perseus, Ophiuchus and Pipe Nebula molecular 
clouds - which all fail at forming massive stars -, and of molecular 
clumps selected for s i gnpos ts of massive star formation activity led 
IKauffmann & Pillai 1 420ich to conclude that the mass-radius space 
of molecular structures is characterized by a threshold 

m(r) = 870M e (r/pc) 1 ' 33 (2) 

for massive star formation (MSF). That is, the mass-size sequence 
of molecular clouds devoid of MSF lies below th e aforemen- 
tioned limiting law (fig. 2a in IKauffmann et al. 2010b), while MSF 
molec ular clumps are located above (fig. 2b in lKauffmann & Pillai I 
l20ld) . In other words, MSF demands a mass of molecular gas en- 
closed within any given projected radius higher than what Eq. [2] 
predicts. In what follows, we refer to Eq. [2]as the MSF l i mit. N ote 
that following the terminology of IKauffmann & Pillai I fcOld) . a 
MSF region comprises HII regions, which have been character- 
ized in literature as regions forming high mass stars, i.e. stars with 
masses > 8A/q. 

The aim of this contribution is to demonstrate that this lim- 
iting law for massive star formation in the mass-radius space can 
be explained in terms of a volume density threshold for overall 
star formation. The body of observational evidence in favour of 
a volume (or number) density threshold for star formation has 
steadily been growing over the past few years. On the observa- 
tional side, there is a tight association between the star forma- 
tion rate and the mass of high density molecular gas (i.e. hy- 
drogen molecule number densities «H2 — 10 4 5 c»r 3 ), in both 
Galactic and extragalactic environments dGao & Solomon 1 12004| ; 
IWu et alj 12001 lLada Lombardi & Alves Il2010|) . On the theoret- 
ical side, lElmegreenl | |2007|) notes that number densities njj2 ^ 
10 5 e/n -3 enhance microscopic effects able to accelerate star for- 
mation significantly (e.g. magnetic diffusion in the molecular gas). 



The effective radius is here defined as the radius of the disc whose surface 
area is identical to that of the contour 



The existence of a number density threshold, for star forma- 
tion also implies that the mean number density of cluster-forming 
regions is a few times n,;,, the exact factor depending on the clump 
density profile. This i s in l ine with the conclusion reached by 
IParmentier & Rroupa I d201lh that constant volume density cluster- 
forming regions are needed to preserve the shape of the young clus- 
ter mass function through the first 50-100 Myr of cluster evolution. 
This effect results from the mass-independent tidal field impact 
experienced by star clusters dynamically responding to the expul- 
sion of their residual star-forming gas when the mean volume den- 
sity of their gaseous precursors is constant. That the shape of the 
young cluster mass function in the present-day Universe remains 
time-invariant for ~ 100 Myr is suggested by many studies (e.g. 
IChandar et alll2010b . [For a more com prehensive summa ry of the 
arguments above, see also section 2 of IParmentier I j20TTb l. To in- 
vestigate how the volume density threshold for overall star forma- 
tio n and the massive star fo rmation limit observationally inferred 
by IKauffmann et alj J2010b1) link to each other is therefore a most 
timely issue. 

The outline of the paper is as follows. In Section [2] we build 
a model that relates the MSF limit to the minimum mass of dense 
star-forming gas needed to form a massive star. Our conclusions 
are presented in Section|3] 

2 MASSIVE STAR FORMATION IN DENSE 
STAR-FORMING GAS 

2.1 Model 

The observed relation between the mass m ec \ of an embedded- 
cluster and the mass m* jnax of its most-massiv e star can be approx- 
imated by e.g. the semi-analytical model of IWeidner & Kroupa I 
fcOOfj. solid thick line in their fig. 1). It shows that the formation of 
stars with masses m* higher than ~ 8-10M,;, requires clusters with 
a stellar mass m ec \ > 100 Mq, Given that a typ ical SFE in cluster- 
forming regions is ~ 30% dLada&Ladal2003l) . the clusters' initial 
gas masses must be in excess of m,i l cr j, ~ 300 Mq for them to form 
massive stars0 Building on the hypothesis that star formation takes 
place in gas denser than a given number density threshold n,/,, to 
understand the origin of the MSF limit therefore equates with un- 
derstanding how the minimal star-forming gas mass, m,/, 
lates to the mass and radius of molecular structures. That is, a 
star-forming gas mass m,/, larger than m,/, cr j t leads to the forma- 
tion of massive stars (/n*. HMX ^ 10A/©), while m f ;, < m t j l cr i, is con- 
ducive to the formation of clusters less massive than m ec i ~ 100M© 
hence to the formation of low- and intermediate-mass stars only 

{nU, m ax < IOMq). 

For the sake of simplicity, let us consider spherically sym- 
metric molecular clumps. The assumption of spherical symmetry 
for molecular structures in general is clearly an oversimplification. 
However, in the case under scrutiny here, we are interested in dense 
molecular clumps hosting forming-star-clusters, rather than the fil- 
amentary giant mol ecular clouds containing them. As an example, 
iBeltran et all d2006l) mapped in 1.2-mm continuum a large sample 
of clumps. They find the mean and median ratios of the full widths 
at half maximum of their clumps along two perpendicular axes x 

2 Note that an SFE of 30% is also of the order of the formation effi- 
ciency required for a star cluster to survive the expulsion of its residual 
star-for ming gas, atbeit largety deple ted of its initial stellar content (see 
fig. 1 in lParmentier & Gilmore |2007|) 
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Figure 1. Comparison between the observed MSF limit and a grid of iso- 
m,i, lines, where m,/, is the mass of star-forming gas (i.e. «H2 ^ nth) con ~ 
tained by a clump of mass r c i lmp and m c i„ mp . Note that these i so-m^ lines 
are of the same nature as those shown in the top panel of fig. 3 in lParmentier I 

feoiih 



andy, FWHM x /FWHMy, to be 1.04 and 0.96, respectively, which 
justifies the assumption of spherical symmetry. 

As for the density profiles of molecula r clumps, power-law s 



have been put forward by various studies, e.g . lHeaton et al. I ll 19931) 
lHatchell et all (bOOOh.lBeuther et al.l fe002l) , lFontani et al. I J2002h 



and Mulleretal. (2002): 



Pclump(s) = kp s p , (3) 

with s the distance from the clump centre and k p a normalizing fac- 
tor. We refer to p as the 'density index'. Studies quoted above find 
1-5 ^ p ;$ 2, In what follows, we will neglect potentially existing 
radi al variations o f p wi thin molecular clumps. 

IParmentierl d201ll) (her eq. 3) shows that the mass of gas 
denser than a volume density threshold, p,/,, contained by a spheri- 
cal clump of mass m c i ump , radius r c [ ump and density index p obeys: 

/ 3 _ p \(3-i0/, 3/ _ 3(3 

th= \4^) m ch,mp r chnnp ■ < 4 ) 

The mass of gas relevant to star formation is m f /,, rather than the 
overall clump mass m c [ ump . In what follows, we refer m,/, as the 
'mass of star-forming gas' (or mcFRg = m th m Eq.QJ. Conceptually, 
molecular clumps in the present paper correspond to picking up a 
single point along a sequence of mass-radius measurements in e.g. 
top panel of fig. 2 inffCauffmann et alj J2010bl) (see also our Fig.©- 

Any convincing model accounting for the MSF limit in the 
mass-radius plane must address two distinct issues: the slope of the 
limiting law and its vertical location (equivalently, its intercept). 
This issue is easily dealt with using Eq.[4) 

Figure Q] shows the lKauffmann et all f2010bl) MSF limit (dot- 
ted black line with asterisks) along with lines of constant m,/,, 
for a density index p = 1.5 and a density threshold for star for- 
mation n t f, = 10 4 cm~ 3 (which equates to p f /, ~ 700M (;i .pc~ 3 ; 
Lada. Lombardi & Alves l l2010l) . It immediately appears that - for 
this sensible choice of parameters - the observed MSF limit agrees 
nicely with the line m,/, ~ 300Mp, at least when r c i ump ^ 0.3 pc. 
mth — 300Mq is precisely the lower bound of the regime for mas- 
sive star formation. At this stage, we can already conclude that the 
hypothesis of a number density threshold for overall star forma- 
tion and the existence of an MSF limit in the mass-radius space 
are clearly linked. Below the MSF limit, star-forming gas masses 
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Figure 2. How the location of iso-m,/, lines in the radius-mass space re- 
sponds to model input parameter variations. Top and bottom panels con- 
sider two distinct masses of star-forming gas, either m,,, = 30Mq or m,^ = 
300.W. . Assuming SFE~ 0.3, the first case corresponds to the formation 
of a IOMq star out of its individual pres-stellar core, while the second case 
is relevant to the formation of a 10M Q star within a cluster assuming the 
re lation between the maximum stellar mass and the embedded-cluster mass 
of lWeidner & Kroupa I f2O06l) . Values for p and are quoted in the key. 



are lower than 300 Mq, which results in clusters less massive than 
~ 100 M©, thus, stars less massive than ~ 1OM0. 

The 'break-point' in each iso-ra,/, line (e.g. at r c i ump = 1 pc 
for in,/, = 5000 M,;,) stems from clumps located to the left of the 
'break-point' having a density at their outer limit at least equal 
to the assumed threshold «,/, (i.e. p(r c i ump ) ^ p t h). As a result, 
m dump = m th f° r those clumps. To the right of the break-points, 
cluster-forming regions (i.e. where «H2 ^ n ih ) represent a fracti on 
only of their host-clump volume (see fig. 2 in IParmentierfeoi ll) . 

Equation|4]explains straightforwardly why the slope of an iso- 
m,i, line agrees well with that of the MSF limit. A constant ra,/, leads 
to m c i„ mp oc r^ lu ^l . For density indices 1.5 ^ p ^ 2, the slope 3 — p 
of an iso-ra,/, line therefore ranges f rom 1.5 to 1, which brac kets 
nicely the slope of 1.33 estimated bvfKauffmann et aT] J2010bl) for 
the MSF limit (see Eq.©. 

The exact location of the inherit line - our model-proxy to 
the observed MSF limit - in the radius-mass space depends on 
our choice of model parameters, density index p and number den- 
sity threshold n,/,. This issue is investigated in Fig. [2] with p = 1.5 
or p = 2, and = 10 4 cm~ 3 or = 10 5 cm -3 . For the sake of 
completeness, we envisage both clustered star formation (bottom 
panel), as detailed above, and individual star formation (top panel). 
As for the latter, we also assume SFE ~ 0.3 based on the compari- 
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so n between the corqj mass function and the stellar IMF performed 
bv lAlves et al .1 120071) for the Pipe dark cloud. As a result, the for- 
mation of an individual massive star requires a mass of star- forming 
gas larger than m,i l cr j, = 30M e . 

The formation of a IOMq star via a clustered mode of star 
formation requires a larger amount of star-forming gas - hence a 
more massive clump - compared to the formation of an individ- 
ual star. As a result, the iso-ra,/, lines move upwards in the bot- 
tom panel with respect to the top one. A higher density threshold 
for star formation implies that the build-up of a given mass m,f, of 
star-forming gas requires a larger amount of molecular gas (since 
a lower mass fraction of the clump takes part in the star formation 
process). As a result, the n,/, = 10 5 cm -3 lines are located above 
their n,/, = 10 4 cm -3 counterparts. Similarly, a shallower clump 
density profile is conducive to a smaller fraction of star-forming 
gas, and p = 1 .5 thus leads to higher clump masses than p = 2. Fig- 
ure[2]demonstrates that all these sensibl e parameter values can ex - 
plain the observed MSF limit inferred by Kauff mann et all l feOlObh . 
possibly with the exceptions of the = 10 4 cm~ 3 , p = 2) model 
for the individual star formation, and the («,/, = 10 s cm -3 , p = 1.5) 
model for the clustered mode of star formation. The difference be- 
tween those models and the MSF limit amounts to almost a factor 
of 10, while in all other cases, the difference is never larger than a 
factor of three. We come back to this point in Section |2~4l 



2.2 Uncertainties on inferring the density threshold 

Conversely to what we have done in the previous section, one could 
infer the density index of molecular clumps and the density thresh- 
old for star formation based on the MSF limit. Since the MSF limit 
equates with a line of constant m,/,, we can match Eq.[2]to: 



'Ichtmp 



(5) 



l j p J clump ' 

which is another form of Eq. |4] This gives p ~ 1.7 and p,/j ~ 
350M:,.pc~ 3 , or ~ 0.5 x 10 4 cm~ 3 . The inferred density in- 
dex p is in excellent agreement with measurements of den- 
sity profile s lopes based on e.g. dust continuum emission (e.g. 
iMuller et alj |2002|) . The inferred density threshold n t f,, how- 
ever, is lower than what is found based on extinction maps and 
young stellar object census in m olecular clouds («,/, ~ 10 4 cm -3 , 
lLada. Lombardi & Alvesll2010l). HCN-based studies of molecular 
gas (nth - 3 x 10 4 cm" 3 , iGao & Solomon Il2004 . and H 13 CO+- 
m apping of molecula r gas (n,/, ~ 10 5 cm~ 3 , bottom panel of fig. 3 
in IParmentier I I20TH) . The cited observations of the high dipole- 
moment species, HCN and H 13 CO + , sample regions of high den- 
sity molecular gas in which the 1-0 transition of these molecules 
is observed. Estimates of the threshold density n,/, for star forma- 
tion are often chosen to be the critical density of these transitions^ 
Therefore, the density threshold for star formation inferred from the 
MSF limit, «,/, ~ 0.5 x 10 4 cm~ 3 , is lower than these line transition 
densities, as well as lower than the density threshold inferred by 
lLada. Lombardi & Alvesl d2010l) . We come back to this point fol- 
lowing a discussion of the various uncertainties affecting the com- 
parison between predicted and observed MSF limits. 

3 We adhere to the following nomenclature: the word 'core' refers to the 
gaseous precursor of an individual star or of a small group of stars, while 
the term 'clump' is designated for regions hosting cluster formation. 

4 The critical density of a spectral line transition is defined as the density 
for which the collisional excitation rate to the transition's upper energy level 
is equal to its radiative decay rate. 
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Figure 3. Mass-size measurements for two clumps (filled red and blue cir- 
cles) with distinct density indices (p = 1.5 orp = 2). Solid symbol-free lines 
depict three-dimensional mass-size measurements, i.e. mir)(r a p). Upper 
dotted symbol-free lines give the two-dimensional mass-size measurements 
m2D(Xap) by accounting for the clump material located in the foreground 
and background of the sphere of radius the aperture r ap , as illustrated in the 
top panel where the line of sight is a vertical. The dotted (black) line with 
open squares describes another type of clump mass-radius relation, namely, 
how clump outer radius and total mass are related. For instance, both clumps 
depicted here have the same mean column density, TV = 2 x 10 22 c nr 2 . This 
type of relation is the main topic of lParmen tier & Kroupa ( 20 Hj) and is not 
to be confused with the radial mass distributions within individual clumps, 
of relevance for the present contribution 



Molecular clump mapping may underestimate the initial gas 
mass since a fraction of the gas mass has been 'fed' to recent 
and ongoing star formation. An SFE~ 0.3 implies that observed 
gas masses m /,, should be increased by a factor ~ 1.5 to re- 
cover the mass of gas 111,1, at the onset of star formation (since 
SFE = m ec i/m th =m ec i/(m ec i +m obs ) and, thus, m th /m obs = (1 - 
SFE)~ l ). However, this correction factor xl.5 applies to the lim- 
ited volume of cluster-forming regions of mass m,/,. The global SFE 
measured over the scale of an entire molecular clump is necessar- 
ily lower (local SFE in cluster-for ming regions versu s global SFE 
in molecular clumps; see fig. 8 in IParmentier II201II) and the cor- 
rection xl.5 for the gas 'lost' to star-formation becomes an up- 
per limit. Stellar-feedback-driven clump gas dispersal constitutes 
another channel through which observed gas masses of molecular 
clumps may be lower than their pre-star formation contents. To as- 
sess the impact of this effect is well beyond the scope of this paper. 
We will assume that, due to star formation and gas dispersal com- 
bined, pre-star formation clump masses may be higher than their 
observed counterparts by at most a factor 2. In other words, ac- 
counting for star formation and gas dispersal implies that a correc- 
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tion factor < 2 should be applied to the intercept of the observed 
MSF limit prior to comparing it to model outputs. However, the 
observed MSF limit is yet affected by another effect - an overesti- 
mating one - of about the same amplitude, as we now explain. 

Observed masses are two-dimensional masses seen through 
an aperture corresponding to a contour of constant surface density 
and, therefore, include material in the foreground and background 
of the spherical region contained within this aperture. This effect is 
depicted in the top panel of Fig.[3]wifh the line-of-sight the vertical. 
The projected mass, m2o{r ap ), seen through the aperture of radius 
r ap obeys: 

= Ankp J s 2 ~ p ds 



i / s 2 r 2 

+ 4**,/ [l 



- p ds, (6) 



where the first and second terms on the rhs account for the sphere 
of radius r ap centered onto the aperture (mjr>(r ap )), and for the 
correction for the background and foreground material (m cor (r ap )), 
respectively. The bottom panel compares the three-dimensional 
mass m iD (r ap ) = m c i ump (r ap /r c i ump ) 3 ~ p (solid lines) to its two- 
dimensional counterpart m2o(r np ) (dashed lines) for two clumps 
of different density indices p. The ratio between both masses de- 
pends on the aperture size compared to r c [ ump , and on the density 
index p. It appears that the observed mass overestimates the three- 
dimensional mass by a factor not exceeding 2.5 (see Kauffmann et 
al., in prep for a more detailed analysis). 

The uncertainties overestimating (projection effects) and un- 
derestimating (star formation and stellar feedback) the MSF limit 
are therefore of similar amplitude, which limits the uncertainties of 
a relative comparison between models and observations, as done at 
the beginning of this section. 

We stress here tha t the mass-radius relations de- 
fined by Kauffmann et al J2010al fhh on the one hand, and 
IParmentier & Kroupa 1 20 111) o n the other hand have distinct 
physical significances. While in iKauffmann et alj feOlOal fhl). the 
mass-radius relation describes the radial distribution of the mass of 
mo lecular gas regions, the mass-r adius relation discussed in depth 
by IParmentier & Kroupa I d20 1 lh in the framework of the early 
survival of star clusters in a tidal field refers to the mass-radius 
relations of populations of clumps, that is, the relation between the 
total mass of individual clumps and their radius at their outer edge. 
For instance, the two clumps of Fig.[3]have identical mean surface 
densities (E ~ 0.1 g.cm~ 2 or — 2 x \0 22 cm~ 2 ) but different 
density profiles (p =1.5 and p = 2). Note that the choice of mass 
and density index for the two clumps in that figure is arbitrary and 
entirely chosen for illustrative purposes. 

When working in the framework of the clustered mode of 
star formation, another uncertainty stems from the maximum stel- 
lar mass vs. embed ded-cluster mass rel a tion. B ased on the semi- 
analytical model o f IWeidner & Kroupa I J200^) , we have assumed 
that m ec i = 100M,: ; leads to m* Jl1tu = 10 M W . Yet, inspec t ion of 
the data points most recently collected by IWeidner et al. I J2010L 
right panel of their fig. 3) indicates that, on the average, IOMq 
stars are hosted by clusters ~ 50M G in mass. This decreases our 
estimate of 

ftith crit down to 150 Mq (again with SFE — 0.3 in the 
cluster-forming region) and increases by a factor 2.4 the inferred 
dens ity threshold, resulting in n,i, ~ 1 4 cm~ 3 , in good agreement 
with lLada. Lombardi & Alves I d2010l) . 
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Figure 4. How to estimate graphically the star-forming gas mass for a 
clump with given mass and radius (large filled (red) circle), m,;, follows 
from the intersection between the clump mass radial distribution and the 
mean number density < ncfn g > matching given number density threshold 
n,i, and density index p (see Eq.|7J. 



Finally, we expect our analysis to break down for large enough 
spatial scales, that is, on scales comprising several clumps since the 
assumption of spherical symmetry then ceases being valid. In the 
Galaxy, our theory thus may not apply beyond radii of several pc. 
Note also that whether the observed MSF limit stands for effective 
radii larger than lOpc, or not, is currently not known. 



2.3 How much star-forming gas? A swift graphic approach 

To estimate the mass m,/, of star-forming gas of a molecular clump 
- hence its ability to form massive stars for a given SFE - can even 
be done graphically, without resorting to Eq. [4] Let us consider a 
clump of mass m c [ ump ~ 3 x 10 4 M Q , radius r c \ lmp ~3pc and den- 
sity index p = 1.5 (large filled (red) circle in Fig. |4)- The solid 
(black) line traces the clump three-dimensional mass radial distri- 
bution, i.e. m(s) oc s^~ p K The mean density of the st ar-forming 
region of mass m,/, obeys (see eq. 5 in Parmentier 1l201lh : 



< n CFRg > = 



-n t h ■ 



(7) 



A density index p = 1.5 is thus conducive to a mean number den- 
sity < ncFRg > within the star-forming region twice as high as on 
its outer bound. This mean density < ncFRg >= 2 x 10 5 cm~ 3 is 
depicted as the (blue) dashed line in Fig. [4] (assuming a thresh- 
old /;,;, = 10 5 cm -3 ). The intersection between the line of constant 
< ncFRg > an d the clump mass radial distribution, indicated by a 
filled (blue) circle, indicates therefore the radius and the mass cor- 
responding to a mean number density < /icf«? >= 2 X 10 5 cm~ 3 
hence a limiting density n,/, = 10 5 c»r 3 . In other words, the mass 
indicated by the small (blue) filled circle defines the mass m,/, 
of star-forming gas, that is, 500 Mq. A reading of the bottom 
panel of Fig. [2] actually indicates that a clump characterized by 
those mass, radius, density index and number density threshold has 
m th ^ 300 Mq . It is thus possible to estimate graphically the amount 
of star-forming gas contained by a molecular clump, by searching 
for the intersection between the clump mass radial distribution and 
the mean volume density of the gas denser than the threshold n,;,. 
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Figure 5. Top panel: Combination of two MSF models: formation of 
a 10M ( :,-star out of a density peak with n,i, = 10 5 cm~ 3 (red dash-dotted 
line), and formation of a 10 M; -star within a CFRg with = 10 4 cm -3 
(blue dashed line). The combined model explains the MSF limit down to 
spatial scales of ~ 0.05 pc, or clump mass ~ 30M e , in excellent agree- 
ment with the smallest clumps se lected for massive star formation activ- 
ity by iKauffmann & PillaTI l201fj) . Note that, for the sake of clarity, the 
pre-stellar core model (red dash-dotted line) has been shifted by -0.1 in 
log(m). Bottom panel: sam e as top panel, but completed with the data of 
IKauffmann & Pillai I <20lTj) and the shaded area discussed in the text. 



2.4 Individual- and clustered-MSF: a combined approach 

Figure [TJ shows that the observed MSF limit and the CFRg model 
in,/, ~ 300 Mq agree reasonably well with each other when r c i ump > 
0.3pc. For smaller clumps, the model corresponds to a line of con- 
stant mass (i.e. all the clump gas is denser than the density thresh- 
old riff, = 10 4 c/n~ 3 ), while the observed limit retains its slope of 
1.3 down to spatial scales of a few 0.01 pc ( see top and middle 
panels of fig. 2 in IKauffmann & Pillai I l201(il) . Top panel of our 
Fig. [2] suggests that the regime r c i, lmp < 0.3 pc can be accounted 
for by the formation of stars out of their individual density peaks 
withn,/, ~ 10 5 cm -3 . This 'individual-star-formation' picture is ac- 
tually expected if observations look into a forming-star-cluster, at 
the spatial scale of individual pre-stellar cores. 

Top panel of Fig. [5] depicts two models, one for the forma- 
tion of a star cluster out of a CFRg of mass m,,, ~ 150Mq and 
density threshold n,/, = 10 4 cm -3 (blue dashed line), the other for 
the formation of a star out of its individual pre-stellar core of 
mass m,i, ~ 30Mq and density threshold n,/, = 10 5 c m~ 3 (red dash- 
dotted line). Assuming SFE ~ 0.3 in both cases dLada & Lada I 
l2003l : lAlves et alj|2007l) . these thresholds in the mass-size space 
of molecular clumps are consistent with the formation of a 10 M Q - 
star on the spatial scale of either a CFRg, or of a pre-stellar core. 



Note that since the SFE enters the m c i ump -r c i ump relation under a 
power p/3 ~ 0.56 (see Eq.[5j, its exact value does not influence the 
inter cept of the models s ensitively. 

Elmegreen ] < l201ll) confirms that both approaches can be 
merged into one single coherent picture. In his model, the average 
density profile of molecular clumps is described as a decreasing 
power-law (his equation 5) which, once convolved with a density 
probability distribution function for supersonically turbulent gas, 
allows to account for the formation of pre-stellar cores by local 
turbulence compression all through a clump. That i s, each molecu- 
lar cl ump contains many pre-stellar cores (see also iMcKee & Tan I 
|2002|) . Pre-stellar cores correspond to density peaks of at least 
n,f, = 10 5 cm~ 3 becaus e stars form fastest where ;iH2 > 10 5 era -3 
(see lElmegreen"! I2007I his section 3.6). The mass fraction of 
clump gas turned into pre-stellar cores by the turbulence is not 
uniform through a clump, but increas es with the clump average 
density (see fig. 1 in lElmegreen 11201 II) . That is, per unit gas mass, 
pre-stellar cores have a greater probability of forming near the 
centre for a centrally concentrated clump. In that respect, our 
assumption of a number density threshold n,,, = 10 4 cm~ 3 for 
cl ustered star f ormat ion constitutes a first-order approximation 
to lElmegreen I ( 1201 ll) 's model, i.e. the increased likelihood of 
forming pre-stellar cores towards molecular clump central regions 
is replaced by a Heaviside Function //(«H2 — 10 4 cra~ 3 ): SFE = 
if n m < 10 4 cm -3 , SFE ^ if n m ^ 10 4 wr 3 . 

Whether all massi ve stars form in star clusters re- 
main s disputed (see e.g. Ide Wit et al. 1 120051 IParker & Goodwin I 



| 2007l: iGvaramadze & Bomans I bOOafschilbach & Roser 1 120081 
IWeidner et al. 1120 lot for contrasting points of view). We note that 
Fig.|5]alone cannot disentangle between the non-clustered and clus- 
tered modes of massive star formation. To illustrate this, let us con- 
sider a clump with r c i, tmp = 1 pc and m c \ ump = 2000 M©. The gas 
masses with «H2 > 10 4 cm" 3 and 11^2 > 10 5 c/n~ 3 it contains are 
800 M,:, and 135 Mq, respectively (using Eq. [4] or the graphic ap- 
proach of Section |2~3l . With SFE ~ 0.3, this leads to the formation 
of either an embedded-cluster of mass ~ 240 M Q , or a single star of 
mass ~ 45 Mq. On the average, the most massive star hosted by a 
240 Mq -cluster has a mass ~ 15 .V/. . Therefore, both models differ 
sensitively by the mass of the most massive star formed within the 
clump. Due to the high density of molecular clumps, most of the 
bolometric luminosity given off by the stars forming within their 
interiors is reprocessed into the infrared. Besides, the luminosity of 
main sequence stars depends strongly on their mass (i.e. L °= m 3 5 ). 
For the case depicted above, the luminosity of the single-star model 
thus outshines that of the cluster model by a factor ~ 40. Com- 
bining our model to measures of molecular clump infrared lumi- 
nosities appear therefore as a promising tool to test whether some 
massive stars can form outside a cluster environment. 

Bottom panel o f Fig. [5] reproduces the data collected by 
IKauffmann & Piliail l20ld) in the middle panel of their fig. 2, 
with the same colour- and symbol-codings. A handful of molec- 
ular clumps are located in the area bounded by the MSF limit and 
the nm, = 30Mp, and m,/, = 150M Q models (left shaded area in bot- 
tom panel of Fig.|5j. Given the small spatial scale (< 0. 1 pc) and the 
high mean density («H2 > 10 5 cm -3 ) of many of these 'clumps', the 
term 'core' is better relevant. Although those cores might indicate 
non-clustered massive star formation, they could also belong to a 
larger - cluster-forming - region with a total mass > 150— 300 Mq, 
hence a region expected to form massive stars. Inspection of the 
spatial surroundings of the cores in the shaded area of Fig. [5] will 
allow to settle this point. 
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3 CONCLUSIONS 
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Building on the model presented in IParmentierl d201ll) . we have 
demonstrated that power-law density profiles for molecular clumps 
combined with a volume density threshold for star formation yields 
mass-size relations for the molecular clumps containing a given 
amount of star-forming gas hence forming star clusters of a given 
mass (Fig. [TJ. We have used this re sult to explain the massive star 
formation limit recently inferred bv lKauffmann & Pillai I d20~ich . 

The mass m c i ump , radius r c \ ump and density profile slope 
p of molecular clumps which contain a mass m t i, of gas 
denser than a given density threshold p,/, obey: >n c [ ump = 

m> tC ( ^i-p ) r clump ( our Eq|5}. That is, this relation pro- 

vides the clump mass-size relation corresponding to any given 
amount m,/, of star-forming gas for a star formation process driven 
by a volume density threshold p,/,. Note that the assumption of 
a density threshold for star formation implies that formation of 
star clusters takes place in a limited region of their host molecular 
clumps (i.e. the cluster-forming region does no t necessarily matc h 
the volume of the whole clump, see fig. 2 in IParmentier I I20TH) . 
Our analysis requires very few further assumptions: (i) molecu- 
lar clumps hosting forming clusters are spherically symmetric and 
(ii) the SFE in cluster-forming regions is mass-independent (as de- 
manded by the mass-independent infant weigh t-loss of young star 
clusters, fig. 1 in [Parmentier & Gilmore 1120071) . 

In the framework of the clustered mode of star formation, we 
link the mass of a cluster to t he mass of its most-massive star with 
the semi-analytical model of IWeidner & Kroupa thereby 
implying that the formation of a ~ lOM,., star requires a cluster 
mass of ~ WO Mr.. . As star formation efficiencies in cluster-forming 
regions are expected to be of order SFE ~ 0.3, the minimum mass 
m th,crit ma t a cluster-forming region must have to form a massive 
star is therefore ~ 300 Mq. 

Armed with this estimate of m f /, cn/ , we have compared the 
corresponding model for the mass-size sequence of molecular 
structures (our Eq. |5) with the observed massive star form ation 
limit m(r) = 870M,7 ; (r/pc) L33 of lKauffmann & Pillai I fj)ich . We 
derive an estimate of p ~ 1.7 for molecular clump density in- 
dices - which is in excellent agreement with other estimates of 
the literature -, and a number density threshold for star formation 
riii, ~ 0.5 x 10 4 cm -3 . Th i s is lo wer than what is inferred by e.g. 
lLada. Lombardi & Alves n t h — 10 cm i ) based on extinc- 

tion maps and young stellar object inventories of molecular clouds. 
While the difference may be related to uncertainties inherent to our 
method, we also speculate that this may indicate that IOMq stars 
can form out of clusters slightly less massive than 100M ; ^, as sug- 
ge sted by the star cluster and stellar masses most recently collected 
bv lWeidneretal.l l l2010l) . 

The formation of massive stars as members of a star cluster 
allows to explain the observed MSF limit down to a spatial scale of 
0.3 pc. For smaller radii of molecular structures, we show that the 
MSF limit is consistent with the formation of a IOM0 star out of 
its individual density peak with n t \ r ~ 10 5 cm . Th is is the density 
of pre-stellar cores predicted by|m megreen 1 (|2007|) . The observed 
MSF limit therefore embodies information about the formation of 
massive stars as star cluster members, and the formation of massive 
stars out of their individual pre-stellar cores. In this framework, the 
density threshold n,h — 10 4 cmT^ for clustered star formation prob- 
ably represents the mean density above which the formation of lo- 



cal density peaks with : 
turbulent gas. 



10 cm is favoured in supersonic ally 
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